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Abstract 

(N . 

We use the strong coupling and hopping parameter expansions to cal- 
culate the pion and rho meson masses for lattice Yang-Mills gauge theories 
with fermions in irreducible two-index representations, namely the adjoint, 
Qh| symmetric and antisymmetric. The results are found to be consistent with 

■ orientifold planar equivalence, and leading order 1/N C corrections are cal- 

culated in the lattice phase. An estimate of the critical bare mass, for which 
the pion is massless, is obtained as a function of the bare coupling. A com- 
parison to data from the two-flavour SU(2) theory with adjoint fermions 
gives evidence for a bulk phase transition at f3 c ~ 2, separating a pure 
lattice phase from a phase smoothly connected to the continuum. 



1 Introduction 



Recently, there has been an ongoing effort to use lattice techniques for 
studying properties of gauge theories beyond QCD, with a large number of 
^ . colours (N c 3> 3) and/or fermions in higher dimensional representations. 

The earliest motivation for this has been the proposed 'orientifold planar 
equivalence' [U [2] which is valid at infinite N c , and provides a link be- 
tween pairs of gauge theories (for a general review, see [3]; for a lattice 
formulation, see [4]). As a particular case, it predicts that 'adjoint QCD', 
'symmetric QCD' and 'antisymmetric QCD' (by which we mean theories 
whose action is the same as QCD but with the fermions in the respective 
representation) all have the same bosonic sector at infinite N c . We refer to 
these theories collectively as 'orientifold theories'. It is interesting that by 
taking fermions in the adjoint, the one-flavour theory is identically A/" = 1 
Super Yang-Mills, and so we can copy analytical predictions obtained using 
supersymmetry to the other two theories py. Furthermore, for N c = 3, the 
antisymmetric theory becomes one-flavour fundamental QCD, suggesting 
a pathway for making real predictions in a close relative to real QCD - 
provided N c = 3 is 'close' to infinity. However, today this question can only 
be addressed by measuring the size of 1 /N c corrections non-perturbatively 
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using lattice methods. In the case of pure Yang-Mills or quenched QCD 
with fundamental fermions, many studies have found that the corrections 
are indeed small (for a review, see [6]). There has been much less work 
on two-index fermions, though studies are now beginning to appear. In 
particular, in [7], a quenched lattice simulation of the quark condensate 
in orientifold theories was carried out, and a comparison with the analytic 
expression from [8] supports the equivalence. Note however that for two- 
index fermions, the quenched theory and the dynamical theory are different 
at infinite N c , so a definitive result has yet to appear. 

Orientifold theories have also gained attention as candidates for Beyond 
the Standard Model physics, as their dynamics is potentially very differ- 
ent to QCD. In particular, there are proposals for 'Technicolor' models 
of Dynamical Electro- Weak Symmetry Breaking where the Higgs is re- 
placed by a composite bound state of strongly coupled higher-dimensional 
fermions. A recent concrete example making use of such a Higgs sector 



is Minimal Walking Technicolor [9], which in our language is just SU(2) 
adjoint QCD with two flavours. Numerical lattice studies have already 
been carried out to determine the non-perturbative dynamics of this the- 
ory nnimnanaimngnginiiiB], and there is m ° untin g evidence 

of a conformal infra-red fixed point, or at least near-conformal behaviour 
- a requirement for the walking scenario. There have also been numeri- 
cal investigations of conformal behaviour in the case of SU(3) symmetric 



Since, however, unlike QCD, there is no experimental data to guide 
the interpretation of numerical results, it is important to learn as much as 
possible by analytical means. A weak-coupling analysis has already been 
performed [25] and, among other conclusions, gives perturbative estimates 
for the ratio of A parameters, A\ at / Aj^, and for the additive renormal- 
isation of the quark mass and of fermion bilinears in this regime. The 
present paper is written in the same spirit, and looks at the opposite side 
of the lattice phase diagram by studying the meson spectrum in the strong 
coupling regime. 

It is clear that a lattice strong coupling expansion has no relevance 
in the continuum limit, however our goal is to establish analytic results 
which, first, will give a starting point for choosing simulation parameters 
and, second, will provide information on the phase structure of the lattice 
theory. In particular, we derive formulae against which numerical data can 
be compared to ensure simulations are not in an 'unphysicak phase (in the 
sense of not having a continuum limit). 

In addition, the meson masses provide explicit observables for which 
to check orientifold planar equivalence at infinite N c . Formally, a general 
proof of planar equivalence that holds to all orders in the strong coupling 
and hopping expansion has already been presented [1], and the results of 
this paper should be considered a special case. The benefit of our direct 
calculation is that it provides explicit expressions for the meson masses at 
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finite N c . 

Section [2] of the paper establishes the notation and describes the strat- 
egy. The calculation is performed in Section O closely following the dia- 
grammatic method of [26]. We enumerate the relevant diagrams and ex- 
plain the generalisation from fundamental to two-index fermions, including 
a description of special cases which fall outside the general analysis. In Sec- 
tion S] the results are discussed and a connection is made with available 
numerical data, finding two distinct regimes consistent with a lattice phase 
and a phase connected to the continuum. The conclusions are summarised 
in Section [5j 

2 Strong coupling and hopping expan- 
sion approximations 

Here we set up the notation and outline the strategy. Discretising using 
Wilson fermions in lattice units (a = 1), the action and Dirac operator is 

S = S g + ^(x)D(x,y)i>(y) (2.1) 
U(x, n)U{x + (i, v)U\x + u, n)U\x, v) + h.c 



x,y 
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"0 x,fi>u 

D(x,y) = 5 xy -K(x,y) 

K(x,y) = 2K^2[p-V(x,f,)5 y , x+ ^ + P+V(x-fi,^5 y>x .^ , (2.2) 

where P^ = (1 ± 7/J/2 are the standard projectors, we have set the Wil- 
son parameter r to 1, and the expansion parameter k is related to the bare 
quark mass by the usual relation 2k = 1/(4 + mo) We do not specify the 
number of flavours beyond stating Nf > 2, as the final result will coincide 
with the quenched result to the order we work to in the hopping expan- 
sion. The gauge part of the action is the standard sum over elementary 
plaquettes, with the links U always transforming in the fundamental rep- 
resentation of the gauge group. In the Dirac operator the links V are in an 
arbitrary representation, which we will take to be either the fundamental 
or a two-index irreducible representation. 

We want to compute the two-point meson correlator in the triplet chan- 
nel, 

G a/ 3y8{x,y) = (ip a (x)^(x)tp' y (y)i/; s (y)) 

= ~ J V[U](detD)D^(x,y)D^(y,x)e s o , (2.3) 

where tj) and tfj' are different fermion flavours and colour indices are con- 
tracted to make colour singlet mesons. There is a well-established proce- 
dure for computing the correlator in the strong coupling/hopping expan- 
sion approximations: 
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1. Expand the quark propagators as a series in k (hopping expansion]^]. 
To each k is associated a factor given by one of the terms of (|2.2p . 
which geometrically corresponds to a link between two adjacent lat- 
tice sites. The full expansion is constructed by summing all possible 
discrete paths linking the two spacetime points x and y. Note that 
backtracking paths are not allowed as we have made the choice r = 1 
for the Wilson parameter in the fermion action, so that P^P^ = 0. 

2. Expand the fermion determinant in powers of k 4 . This inserts an 
arbitrary number of closed four dink plaquettes (the product and trace 
of four factors (|2,2p ) starting at any point on the lattice, in any 
orientation. There are k (possibly overlapping) plaquettes at order 
k. 

3. Expand the gauge part of the action in powers of 1/<?q (strong cou- 
pling expansion). This also produces four- link plaquettes as in the 
previous step, but containing only fundamental gauge links U and no 
spin factors. 

4. In the total expansion there will then be a mixture of links in differing 
representations (U in the fundamental from the gauge part and V in 
a two-index representation from the fermion part). The only non-zero 
contributions are those where the gauge links come in combinations 
which survive the Haar integrals. We will only need the two-link 
U(N C ) integral, 

j V[U]R a [U] ij R b [U% = ^ ab S u S jk , (2.4) 

where R a [U] denotes the link U in representation a, and dji is the 
dimension of the representation. The difference between U(iV c ) and 
SU(-/V C ) will be discussed in Section [331 

To make progress, we restrict ourselves to a fixed order in k and 1/<7q. 
As a starting point, expand the gauge action e Sg to lowest (zeroth) order 
in 1/<7q, where there are no gauge plaquettes, and all gauge links have to 
come from the fermionic part of the action. This means that in order to 
satisfy (|2.4p . the two paths from two quark propagators must be colinear 
at all points (and in opposite directions jl, see Figure HJ The full path can 
be built up recursively from an elementary building block, M, which to 
this order is just two colinear links. We can then go to higher order by 
generalising M to include non-colinear segments (i.e. to the next order, 
gauge squares, see Section l3~2l) . This method for constructing all paths 
was used in [26J to compute the pion and rho masses in standard QCD. 
We refer to this paper for the technical details, and here only establish the 
notation. 

1 Note that this is a large mass expansion, immediately excluding using this approach for 
studying the open question of whether a theory is conformal in the IR. 

2 There is also a possibility of non-colinear paths, with (|2.4I) satisfied by insertions of pla- 
quettes from the fermion determinant, but this is higher order in k than we consider here. 
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Figure 1: In order for Haar integrals to be non-zero, links must be oppositely 
paired between each lattice site. If we set e~ Sa ~ 1 and detD ~ 1, the only 
possibility is for the quark and antiquark propagators to be colinear. 

Define Gi jCf ^ 7 ,5(n, 0) to be the sum of all paths of length L from the 
origin to n. The full correlator is then 

oo 

G a/ 3 7(5 (n,0) = ^G LiCt/ 3 75 (n,0) . (2.5) 
L=o 

To calculate it, we can relate Gl to Gl~i through 

0) = J2 M ^r (n, n')GL-i^5 T f3(n' , 0) , (2.6) 
n',<r,T 

where M contains the factors describing the propagation between sites n 
and n'. It can then be shown (e.g. by looking at (|2.6[) in Fourier space) 
that the full propagator is then just the sum of an infinite geometric series, 

oo 

G aj a 7 5(n,0) = G L<a i3-ys(n) = (5 a s5-y/3 - M a/37(5 (n, 0)) _1 . (2.7) 
L=o 

We can look at the desired channels by inserting the appropriate T matri- 
ces, 

G A B{n, 0)=Y, (^A)ap(r B ) lS G a p 7 s(n, 0) , (2.8) 

with the r for each of the scalar, pseudoscalar, vector, axial vector and 
tensor channels given by 

rs = - t p = Its r Vfi = -t m 

i 1 

T A„ = -7m75 T Tltu = 7~7|>^ - IvHv , (2-9) 

and normalised as 

TrT^rfi = 5aB = flj h^vjSfjtv, -ji^an^Pv — Sau8f3n)) • (2.10) 

The masses are then extracted by looking in momentum space for poles 
in Gab (or, equivalently, zeros in its inverse) for a particle at rest with 
momentum = (im,0). 
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3 Computing the masses 
3.1 Leading order 

In the framework of section O the problem reduces to writing an expression 
for M(n,n'). For the lowest order colinear paths, the expression is simple: 



M a ^ 5 (n,n>) = (2 K ) 2 ^[(P M ) a5 (P+)^., 



n+fj,,n' 



+ {P^)as{P^8 n -„,n'] • (3.1) 

We have omitted the links V, as the colour contribution does not depend on 
the shape of the path and can be factorised and easily treated separately, 
giving a constant factor dp. Substituting (13. ip into (12. 7|) and taking the 
Fourier transform, we obtain the expression 

r d 4 p ■ - 
G af 3 7 s(n,0) = -d R J j—^e^Gafasip) 

+(P^)as(Pp^e-^} . (3.2) 

Contracting (|3.2p with T matrices then gives the momentum space corre- 
lator. Detailed expressions can be found in [26]; we simply quote the final 
result after substituting p^ = (im, 0), for the pseudoscalar channel Tp, Ta 
and the vector channel Ty i , Tj- 0i (where i is a spatial index) . 

/ 1 — 12k 2 — 4K 2 coshm 7r 4iK 2 sinhm 7r \ , . 

PA ° \ — 4?K 2 sinhm 7r 1 — 4K 2 coshm 7r J 



The zeros of these matrices, and thus the masses, can easily be found by 
solving det(5p^ o = 0, detGy.j, , = 0, to give the well known results 

(1 - 16k 2 )(1 -4k 2 ) . x 

coshm, = l+ y >\ ; 3.5 

(1 - 12k 2 )(1 - 8k 2 ) , . 

coshm p = 1 + —rj- — ^ 3.6) 

8k z (1 — ok z ) 

where, in analogy to QCD, we call the pseudoscalar particle it and the 
vector p. To this order, there is no dependence on representation. 

3.2 Next to leading order 

To write an expression for M(n, n') at higher order one must jointly con- 
sider the l/<?o an d k expansion. We will choose to truncate at order k 6 , 



, 1 — 8k 2 — 4/-t 2 coshm„ 4i\ / i/i 2 sinhm , 

Gy} Tai = \ r P V2 P |(3.4) 

— 4iW ^K 2 sinhmp ^(1 — 4k 2 — 4K 2 coshmp 
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while the order for the strong coupling expansion will depend on the rep- 
resentation. 

We consider the fundamental representation first. This case was al- 
ready studied in [26J, but it will be useful to present it to explain how to 
generalise to higher representations. When expanding e Sg to first order, 
we are allowing the placement of one gauge plaquette on the lattice. The 
edges must be paired with oppositely oriented links, and the only possible 
insertion is within a square of links from the fermion lines. M(n, n 1 ) will 
therefore contain terms propagating from one site to the next not only 
colinearly but also leaving at most one empty square, to be filled by a 
plaquette from the strong coupling expansion. It is also possible to fill the 
square with a plaquette from the fermion determinant, but this introduces 
an extra factor of k 4 , which, together with at least four other fermion lines, 
would be higher order than we consider. There are a number of ways and 
orientations to construct this square, and they are all listed along with their 
order in Table HJ The colour integral for each diagram can be evaluated, 
and all diagrams give a factor l/g^ for the fundamental representation. 

What happens if we change the representation of the fermions? The 
first thing to notice is that the diagrams as described above are all zero 
because of the orthogonality condition (|2.4p . since the fermions are in a 
two-index representation and the gauge links are in the fundamental. Since 
there are no other possibilities at order l/g$, this whole order is zero and 
we must consider the next one. At second order, we can place two gauge 
plaquettes on the lattice. Overlapping them inside the fermion square 
(in the correct orientation, which is representation dependent) will yield a 
non-zero contribution. 

For concreteness, consider the colour factor for the symmetric represen- 
tation. The gauge integrals can be done by writing the reducible product 
of two overlapping fundamental links in terms of symmetric links S and 
antisymmetric links A, 

U ac Ubd = S(ab),(cd) + A [ab],[cd] ( 3 -7) 

where (ab) indicates the symmetric combination, with b > a, and [ab] 
indicates the antisymmetric combination, with b > a. Defining a new basis 
such that S^^^) = SijE l ab E J cd (and similarly for the antisymmetric) with 
i,j running from 1 to da, the gauge integrals take the form 

J-DmSu + A^St (3.8) 

By orthogonality, the second term is zero, and the first is given by (|2.4p . 
The argument is identical for antisymmetric fermion links, with the first 
term being zero and the second given by (12. 4p . Thus, we are back in 
a case exactly analogous to the fundamental, and the calculation follows 
through identically but with a different prefactor owing to the higher order 
expansion of e Sa . There is a factor 2! from the Taylor expansion, and a 
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Spin contribution, including colour prefactor 


1/90 


K 





2 


— < — 
— > — 


(2^) 2 (P-M^ 


Q 


4 




> 











1 


6 


/ 


< 

> 

> 

\ 

< 


+ A 

V 




e| 1 

> 






Q 


6 




IS K P m )aU (( Wft + ( W^" WW] 


Q 


6 


c 




S [( W + {p;p,KPtP;U) w&i 


1 


4 


A 




V 





Table 1: Diagrams contributing to M(n,n'). The first five rows represent diagrams linking sites n to n + /t; there are 
analogous diagrams linking n to n — (i. The last row links n to n + ft + z>; there are also all permutations of terms with 
fi <-> — /t and z> <-> — z>. In the fundamental representation q = 2, while for two-index representations q = A. 



combinatorial factor 1 as there is exactly one way to place the two gauge 
plaquettes. This gives an overall colour factor of l/(2<7 4 ) f° r the symmetric 
and antisymmetric representations. 

The case of the adjoint is similar, except for the fact that the two 
inserted gauge plaquettes must run in opposite orientations. There are 
now two ways to do this, so there is an extra factor of 2 which cancels the 
2! from the Taylor expansion. The result for the adjoint is thus l/g$- 

Having obtained the colour factor, the calculation for the spin part is 
the same for any representation, namely the same as the calculation in 
|26j . Constructing M from the diagrams in Table Q] and extracting the 
poles gives, for the pseudoscalar channel, 

G p \ = ( 5 PP 5 PAo ^ (3.9) 
PA ° V H A P H AoAo J K ! 

H PP = l-4K 2 [3 + 12e fi K 2 -96e/jK 4 

+(1 + 6e R K 2 — 24e R n )coshm 7r ] 

HpA = —Ha p = 4iK 2 (l + 6e R K 2 — 24eRK 4 )sinhm 7r 

H AoAo = 1 + 4k 2 [3 + 9e R K 2 - (1 + 6e R K 2 + 24e jR K 4 )coshm 7r ] , 

and for the vector channel, 

Gy\ = ( ? v * y * ^ Toi ) (3.10) 

H ViV . = ^•[l-4 K 2 (2 + 6e i jK 2 -48e^ 4 

+(1 + 6e R K 2 - 24eRK 4 )coshm p )] 

H Vi T oj = H Toi Vi = -4i8ij-K 2 (l + 6e R K 2 - 24e jR K 4 )sinhm p 

H ToiTi0 = i[l-4 K 2 + 8e i? K 4 + 96e i? K 6 

— 4k 2 (1 + 6e R K 2 — 24e R K )coshm p ] , 

where the constant e R depends on the representation in the following way: 

1 

l 

1 

2d R g^ ' 
1 

Solving detGp^ = 0, detG^. = gives the masses, 

1 - 20k 2 + 64k 4 - 48e fl K 4 (l - 8k 2 + 64k 4 ) 



fundamental e R = ^ (3-H) 

" To 

adjoint e R = - — j, d R = N 2 — 1 (3.12) 

% 

symmetric e R = ^j— j, d R = N C (N C + l)/2 (3.13) 

?Q 

antisymmetric e R = — — j, d R = N C (N C — 1)/2 (3.14) 



cosh m w = 1 + 



8k 2 (1 - 6k 2 + 6e R K 2 (l - 10k 2 + 48k 4 )) 

l - 20k 2 + 96k 4 - 12e fl K 4 (5 - 84k 2 + 384k 4 

coshm = H — rf- — — (3.15) 

8k 2 (1 - 6k 2 + 6e fl K 2 (l - 11k 2 + 48k 4 )) 
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To this order, the meson correlators in the other channels do not contain 
poles leading to real- valued masses. 



3.3 U(N C ) vs SU(N C ) 

The masses as calculated in this section only make use of the U(iV c ) inte- 
gral (|2.4p . so the results strictly apply only to U(iV c ) gauge theories, not 
SU(iV c ). If one is interested in the large- N c limit, this is not a problem 
as the singlet part of U(iV c ) = SU(iV c ) x U(l) decouples, and both groups 
give the same result. However, for small gauge groups, one must also take 
into account contributions of the form 

J 'D[U]Ui 1 j 1 Ui 2 j 2 ■ ■ ■ Ui Nc j Nc = j^^i 1 i 2 ---i Nc € jij2---jN c ■ (3.16) 

With all links expressed in terms of the fundamental representation, the 
diagrams enumerated in Table [T] allow for up to four superimposed gauge 
links - two from the two-index fermions and another two from the insertions 
of (up to) two plaquettes. Equation (|3.16|) is non-zero for N c superimposed 
links, so this integral will contribute for N c = 2,3,4. For N c > 5, to this 
order in the strong coupling expansion, SU(-/V C ) coincides with U(iV c ). In 
addition, representations for some small gauge groups are equivalent to 
each other, and need to be considered as special cases. We look at each 
individually: 

• SU(-/V C ) adjoint — For the adjoint representation, the full contribu- 
tion is captured by including the two orientations of the gauge pla- 
quettes, as we have done in Section [321 and equation (|3.16p plays no 
role. The result (|3.12p is therefore unchanged for all N c . 

• SU(2) antisymmteric — This representation is just the singlet, so the 
theory is simply the free fermion theory. A strong coupling expansion 
is meaningless in this case, and we discard it completely. 

• SU(2) symmetric — For SU(2), the symmetric and adjoint represen- 
tations are unitarily equivalent. As the result for the adjoint has 
already been argued to be correct, the result for the symmetric must 
be the same. Alternatively, we can work directly in the symmetric 



representation, adding the two diagrams in Figure 2(a) The two di- 
agrams turn out to be equal, and this provides the factor of 2 needed 
to give e R = l/(d R gfj). 

SU(3) antisymmetric — This representation is unitarily equivalent 
to SU(3) fundamental. The correct result must therefore be (|3.1ip . 
namely e R = 1/(-/V c <7q). The difference comes about because the 
orthogonality condition (12. 4p does not vanish to order 1 as it does 
for the other two-index representations. 

SU(3) symmetric — From (]3.16p . the only extra three-link diagram 
which could contribute is shown in Figure |2(b)| However, in this 
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+ 



(a) Non-zero contributions to SU(2) symmetric and 
SU(4) antisymmetric. The diagrams are also present 
in SU(4) symmetric, but the one on the right van- 
ishes. 



(b) This diagram for 
SU(3) symmetric van- 
ishes due to symmetry 
properties of the repre- 
sentation. 



Figure 2: Extra diagrams appearing in SU(iV c ), not present in U(iV c ). The solid 
lines are coming from the fermions, and are symmetrised or antisymmetrised. 
The dotted lines are insertions of plaquettes in the fundamental representation. 
All the diagrams which include squares in Table [1] have analogues to the above. 



representation, the symmetrisation of the indices leads to the vanish- 
ing of the diagram. Thus (|3. 13|) is correct without modification for 
SU(3). 

SU(4) antisymmetric — This case is similar to SU(2) symmetric: 
the same two diagrams in Figure 2(a) contribute and they are both 
equal, giving an extra factor of 2 compared to U(4). The result is 

SU(4) symmetric — As for SU(3) symmetric, the symmetrisation of 
the indices make the diagram on the right in Figure 2(a) vanish, so 
(|3.13p is valid as it stands. 



4 Discussion 

The equations (|3.15p give analytical predictions for the pion and rho masses 
in the strong coupling and hopping parameter expansions. We can use 
them to make a few general observations which may help in future lattice 
studies of orientifold theories. The critical value of k where the pion van- 
ishes can be calculated as a function of the bare coupling and is found to 
be 

*> 4 i 1 - 1«) ■ <"> 

Thus, moving away from the infinite coupling limit has the effect of reduc- 
ing k c below 1/4 (although the first order correction is very small). Note 
also that at k c the pion mass is zero but the rho mass remains finite. In- 
deed, the rho mass is always above the pion mass (FigureED, and the strong 
coupling phase is qualitatively similar to fundamental QCD. Furthermore, 
at any given k, there is a definite ordering in the masses (both for the pion 
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-1.5 -1 

Bare quark mass, m. 



■0.5 



Figure 3: Pion and rho mass as a function of the bare quark mass (adjoint 
representation). Parameters: N c = 6, (3 = 10.0, with beta defined in the usual 
way, (3 = 2N c /gQ. The choice N c = 6 was chosen to keep away from the special 
cases discussed in Section [3731 

and the rho) ; the different factors €r are such that the symmetric is always 
heaviest, followed by the adjoint and then the antisymmetric (Figure HJ). 
While this has no physical significance, it is relevant for numerical simu- 
lations, as it tells us that the values of k needed to approach the chiral 
limit will be similarly ordered, with symmetric highest and antisymmetric 
lowest. 

It is interesting to compare (|3.15|) with numerical lattice data. Figure 
[5] plots the pion mass for the SU(2) adjoint theory, using data from [12] 
for two-flavour dynamical simulations supplemented by our own quenched 
simulations. The strong coupling line is plotted for j3 = 0.5, but in practice 
does not move significantly in the range f3 = 0.5 — 3.0. For (3 = 0.5, 
deep in the strong coupling phase, the lattice data falls on top of the 
strong coupling prediction (note that this is not a fit as (I3.15j) has no 
free parameters). Increasing to j3 = 1.5, still quenched, leads to a small 
deviation, more notable for small masses, which is likely to be explained 
by going to higher order in the hopping expansion. At large masses, the 
weak dependence on go seen in (|4.ip is borne out in this phase. The 
dynamical simulations are slightly puzzling: for /3 = 1.5, one would expect 
the quenched and dynamical results to coincide for large masses, as the 
effects of the fermion determinant become negligible. This is not observed, 
suggesting either that the mass is simply not large enough, or that there 
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Bare quark mass, m 



Figure 4: Pion mass as a function of the bare quark mass for the three irreducible 
two-index representations. Parameters: N c = 6, (3 = 5.0. 

could be a small underestimated systematic error in the data. With this 
uncertainty in mind, we can say that the data for 0.5 < f3 < 1.75 is good 
agreement with the strong coupling prediction. In contrast, for (3 > 2, there 
are significant departures from strong coupling, both in the magnitude of 
the masses and even in the qualitative behaviour as one approaches light 
quark masses. This is consistent with the finding in [12J that there is a bulk 
phase transition at f3 c ~ 2, with a strong coupling lattice phase possessing 
no continuum limit for (3 < j3 c , and a phase smoothly connected to the 
continuum for (3 > (3 C . 

The results are also consistent with orientifold planar equivalence, as 
the factors 6r for the two-index representations all tend to the same value 
as N c — > oo (notice that there is a cancellation in 6r between factors of 2 
of different origin, coming together to ensure the equivalence works). In 
the large- N c limit, we have the asymptotic forms of (|3.15p . 





Note, however, that this leading correction to the large-iV c limit is not 
expected to be universal, and may be specific to the regime of validity of 
the expansion: the lattice strong coupling and large mass phase. 
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Figure 5: Comparison of numerical lattice data with strong coupling prediction. 
The quenched data is our own, while the dynamical has been taken from [12J. 
The strong coupling curve is plotted for (3 = 0.5, but in practice shifts very little 
in the range j3 = 0.5 — 3.0. 

5 Conclusions 

We have computed analytic expressions for masses of the tt and the p 
mesons in the strong coupling and hopping expansion (large mass) approx- 
imations, for fermions in the three irreducible two- index representations. 
In the limit N c — > do, the three converge to the same value, as predicted by 
the formal proof of orientifold planar equivalence on the lattice presented in 
[3]. In addition, we have extracted the leading 1/N C corrections, which in 
the strong coupling phase are expressed only in terms of the dimensionality 
of the representation, dn (and the bare quark mass). 

The results are already useful in understanding the lattice phase struc- 
ture emerging from Monte-Carlo simulations. By comparing with recent 
numerical determinations of meson masses in two-flavour SU(2) adjoint 
QCD, we find evidence of two phases, supporting the conclusions of [12]. 
The theory has a bulk phase transition with a strong coupling lattice phase, 
having no continuum limit, on one side (/3 < 2), and a phase smoothly con- 
nected to the continuum on the other {(3 > 2). As more simulations are 
performed for as yet unstudied theories, it is our hope that the results of 
this paper will help in recognising the phase structure and locating the 
correct region of parameter space to use for extracting continuum physics. 
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